FERMILAB-PUB-94 /074-T 

NUHEP-TH-94-6 
March 1994 



J ftp Production from 
Electromagnetic Fragmentation 
in Z° decay 

Sean Fleming 
Department of Physics and Astronomy 
Northwestern University, Evanston, IL 60208 
and 

Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510 

Abstract 

The rate for Z° — ► J/ifi + is suprisingly large with about one event for every 

million Z° decays. The reason for this is that there is a fragmentation contribution 
that is not suppressed by a factor of M^/M|. In the fragmentation limit Mz — ► oo 

with E^p/Mz fixed, the differential decay rate for Z° — > J/ip + £ + £~ factors into 
electromagnetic decay rates and universal fragmentation functions. The fragmentation 
functions for lepton fragmentation and photon fragmentation into J ftp are calculated 
to lowest order in a. The fragmentation approximation to the rate is shown to match 
the full calculation for greater than about 3Mw,. 



Introduction 

Fragmentation is the decay of a high transverse momentum parton into a collinear 
hadron. The differential cross section for the inclusive production of such a hadron in e + e~ 
annihilation factors into differential cross sections da for the production of large transverse 
momentum partons and fragmentation functions D(z) Jl|. The fragmentation function gives 
the probability for the splitting of a parton into the hadron with momentum fraction z. 
These functions are independent of the subprocess that creates the fragmenting particle, 
and can be evolved to any scale via the Altarelli-Parisi evolution equations. 

It has recently been shown that it is possible to calculate the fragmentation functions 
for heavy quarkonium states using perturbative quantum chromodynamics (QCD) 0. Frag- 
mentation functions for several of these states have been calculated explicitly [fj, || £|, [5], ||. 
In particular the paper of Braaten, Chueng and Yuan on charm quark fragmentation || is 

most relevant to the work presented here. Their analysis focuses on decays of the Z° into 
hadrons showing that charmonium production is dominated by charm quark fragmentation. 
In an analagous manner it is shown here that decays of the Z° into charmonium plus electro- 
magnetic particles (leptons and photons) are dominated by electromagnetic fragmentation. 

The process Z° — > ^ + £ + £~ , which is of order a 2 , where a is the electromagnetic 
coupling constant, has a branching ratio of 7.5 x 10 -7 . This is an order of magnitude larger 
than the order-a process Z° — > ipq ]/]], which has a branching ratio of 5.2 x 10~ 8 . That 
makes Z° — > ip + £ + £~ the dominant ip production mechanism in electromagnetic Z° decays. 
The unexpectedly large rate can be explained by a fragmentation contribution which is not 
suppressed by a factor of Mj/M§ |J. In this paper, it is shown that in the fragmentation 

limit Mz — > 00 with E^/Mz fixed, the rate for Z° — > ip + £ + £~ factors into subprocess 
rates for electromagnetic decays and individual fragmentation functions . At leading order 
in a there is a contribution from the fragmentation function for a lepton to split into 

a tp, and a contribution from the fragmentation function for a photon to split into a 

ip. The fragmentation functions are calculated in a manner that is independent of the hard 
process that produces the fragmenting parton, and the fragmentation calculation is shown 
to match the full calculation for greater than about 3M^. 
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The decay rate for Z Q -> tjj + £ + £~ 

The decay rate for Z° — > ^ + was calculated in a model-independent way 

in Ref. |8], |J using the Feynman diagrams in figure la. The diagrams in figure lb also 
contribute to this process at the same order in the electromagnetic coupling a. Fortunately 
they can be neglected. 

In order to understand why only the diagrams of figure la need to be considered it is 
necessary to understand why the rate for Z° — > ^7 is smaller than the rate for Z° — > ip+£ + £~ . 
Naively one would expect T(Z° — > ipj) to be larger than T(Z° — > ip + £ + £~) since the latter 
rate is suppressed by a power of a compared to the former rate. However upon closer 
examination this expectation turns out to be untrue. To see why compare the diagrams 
that contribute to the decay rates T(Z° — > ^ + £ + £~) figure 1 and T(Z° — > ^7) figure 2. 
It is important to note that the lepton propagator in the diagrams of figure 2 is always of 
order 1/M§. In contrast there is a substantial region of phase space where both the photon 
propagator and the lepton propagator in the diagrams of figure la are of order 1/Ml. Thus 
these diagrams will be enhanced by a factor of M|/Ml, compared to the diagrams of figure 
2. The factor of M|/M| is large enough to overwhelm the extra power of a making the 
diagrams of figure la more important than the diagrams of figure 2. The lepton propagator 
of the diagrams in figure lb is always of order 1/M§, so these diagrams are suppressed by a 
factor of a compared to the diagrams in figure 2, and may be neglected. 

The remaining diagrams that contribute to T(Z° — > ip + £ + £~) at the same order in 
a are obtained from the diagrams of figure 1 by replacing the photon propagator with a 
Z-boson. Then the boson propagator in the diagrams of figure la is always of order 1/M§. 
Similarly the lepton propagator in the diagrams of figure lb is also always of order 1/M|. 
Thus these diagrams are suppressed by a factor of a compared to the diagrams of figure 2 
and can be neglected. 

Keeping only the diagrams in figure la, and neglecting the lepton mass, the result of 
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the calculation of T(Z° -> + £+£') is 



r = 4a 2 4 r(z° -> t, 

1 /(y-2) 2 + ^ 



y 2 — w 2 



1+A ry/y 2 -4\ 

dy / g(y,w,X) 

2^A J-^/y 2 -iX 



(y 2 — u> 2 ) s 



y 2 — w 2 



(1) 



where A = M|/M|, y = 2P ■ Z/M 2 Z and iu = 2(p+ - p_) ■ Z/M§. Here P, p_, p + , and Z 

are the 4-momenta of the ip, £~, £ + , and Z°. The parameter can be determined from 
the electronic width T e + e - of the ip to be 
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Using r e + e - = 5.4 keV the photon-to-?/> coupling is g 2 = 0.008. Integrating the function 
g(y,w, A) over w yields the full differential decay rate 

dT 



dE nl 
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where yi = y/y 2 — 4A. In the center of mass frame y = 2E^/Mz- Later on these results will 
be compared to the fragmentation calculation in the fragmentation limit Mz — > 00 with 
E^p/Mz fixed. In this limit Eq. (H) reduces to 



dT „ „ T(Z° 
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The Fragmentation Contribution to Z Q Decay 

In reference |3[ the general form of the fragmentation contribution for the production 
of a ip of energy in Z° decays is given as 



dT(Z° -> ip(E^) +X)=J2jdz df(Z° -> i(^/z) + X, /i 2 ) /i 2 



(5) 



where the sum is over partons and z is the longitudinal momentum fraction of the ip relative 
to the fragmenting parton. All of the dependence on the if) energy has been factored into 

the subprocess decay rate T and all of the dependence on the ip mass has been factored 
into the fragmentation function D^^z, /i 2 ). A factorization scale [i has to be introduced 
in order to maintain this factored form in all orders of perturbation theory. This general 
form was developed in the context of QCD, but it applies equally to QED, where the only 
partons are leptons and photons. This simplifies the general electromagnetic fragmentation 
contribution to 

2 / dz dE £ —(Z°^£~(E e )+X^ 2 )D £ ^(z^ 2 )5(E^-zE e ) 
Jo J dEi 

rl r AT 

+ J o dzj dE^ — (Z° -> 7 (£ 7 ) + X,[i 2 ) D^(z,fx 2 ) S(E^ - zE n \ (6) 

where X are electromagnetic final states, and Eg and E 1 are the lepton and photon energies. 
The factor of 2 accounts for the fragmentation contribution from both the i~ and the £ + . 
Large logarithms of E^/n in the subprocess decay rate T can be avoided by choosing \i on 
the order of E^. The large logarithms of order E^/M^ which then appear in the fragmenta- 
tion functions can be summed up by solving the Altarelli-Parisi evolution equation. In the 
electromagnetic case the evolution of the fragmentation functions is of order a and may be 
neglected. 

It is easy to count the order of a for the leading order fragmentation contributions 
to Eq. (0). The subprocess rate for Z Q — > $r + X is of order 1, while the subprocess rate 
for Z° —>■ 7 + X is of order a. The fragmentation function for a lepton to split into a ip will 
be shown to be of order a 2 , while the fragmentation function for a photon to split into a ip 
will be shown to be of order a. Therefore both fragmentation processes will contribute to 
Eq. (|^) at leading order in a. 

At lowest order in a it is possible to simplify things. The energy distribution for the 
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subprocess Z° — > £ (Eg) + X at lowest order is 

,/! f(Z°^+r)<5(^-^). (7) 



dE t 

Furthermore the photon fragmentation function at lowest order can be written as 

D^z,/! 2 ) = P 7 ^ S(z-l) (8) 

where the function P 7 _>,/, is the probability for a photon to split into a tp, and z is the 
longitudinal momentum fraction of the ip relative to the 7. These simplifications reduce the 
fragmentation contribution to the energy distribution Eq. (|6|) at leading order in a to: 

,/J 'z° -> ^) + = 



^ r(z°-,£ + r) + ^-(z°^ 7 (^) + £ + r, /U 2 )p 7 ^. (9) 

The physical interpretation of the first term on the right hand side of Eq. (H) is that the Z° 
decays into two leptons, with energies Mz/2 on a distance scale of order 1/Mz- Subsequently 
one of the leptons decays into a collinear lepton and ^ona distance scale of order l/M^. 
The physical interpretation of the second term is that the Z° decays into two leptons and 
a photon with energy on a distance scale of order 1/Mz, and the photon fragments into 
afona distance scale of order Note that at this order the only dependence on the 

factorization scale is in D^^p and in the subprocess rate for Z° — > 7 + £ + £~ . 

Given the general form of the fragmentation contribution at lowest order in a in Eq. 
(H), it is only necessary to calculate the fragmentation function D(_^(z, (i 2 ), the fragmenta- 
tion probability P 7 _^, and the subprocess rate for the Z° — ► 7 + £ + £~ . 

Photon Fragmentation 

The fragmentation function D '^^(z , fi 2 ) for a photon to split into a if> can be calcu- 
lated in a manner that is independent of the process that produces the fragmenting photon. 
The Feynman diagram in figure 3a represent such a process at lowest order in a. An un- 
known vertex, represented by the circle, radiates a photon which fragments into a ip. The 



fragmentation probability P 7 _>,/, can be isolated by dividing the cross section ai, for the pro- 
duction of a ip with energy by the cross section cxo, for the production of a real photon 
of energy E 1 = E^, in the limit E^ 3> M$ where fragmentation dominates. 

The general form of the photon production cross section o" is 

^ = FIu~X~ l ^ dPout ^ ( 27F ) 4(54 fe« - k ~ Pont) E l^0| 2 (10) 

where pi n is the sum of incoming 4-momenta, k is the photon 4-momentum, and p ou t is 
the sum of the remaining outgoing 4-momenta. Here [dk] = d s k / (16ir 3 ko) is the Lorentz- 
invariant phase space for the photon and [rfpout] is the Lorentz invariant phase space for the 
remaining outgoing particles, Flux denotes the incoming particle flux for which no explicit 
expression is needed since it will cancel the same factor when <j\ is divided by a . The 
amplitude A for the process can be calculated from the Feynman diagram in figure 3b 

A) = e M , (li) 

where T M is a vertex factor for the production of the photon, for which the explicit form is 
not needed. Squaring and summing over final spins gives 

Y,\A \ 2 = -m;. (12) 

The general form of the 7 — > ip cross section a% is 

vi = j^JldP][dp out }(2Tr) 4 6 4 (p m -P- Pout ) J2\ A i\ 2 ( 13 ) 

where P is the ip 4-momentum. Flux, p in , and p out are the same as described after Eq. fllPP- 
The amplitude A\ can be calculated from the Feynman diagram in figure 3a 

M = 9 ~f ^ {-g^) tr[^], (14) 

where e is the ip polarization vector, and is the vertex factor for the production of a 
virtual photon with invariant mass M^. Squaring the amplitude and summing over final 
spins gives 

El^l 2 = -4vr4ar^r; (15) 
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The vertex factor in Eq. ( |T^D only differs from the vertex factor in Eq. (|TT| ) in one respect, 
it is evaluated at the point k 2 = M 2 instead of k 2 = 0. Aside from the vertex factor the 
tp mass enters the cross section <j\ in the phase space integral. The contribution to photon 
fragmentation will come from a region of phase space where ^> so that the ip mass 
can be neglected. Then the cross section o\ can be written as 

o x = 4:iig 2 a a . (16) 

The fragmentation probability P 7 _>^ can be read off from Eq. (|T6|) : 

P^ = Ang^a. (17) 

Evaluating this numerically gives a fragmentation probability of P^-,^ = 7 x 1CT 4 . 

Photon Subprocess 

The energy distribution of the subprocess Z° — ► 7 + £ + £~ is calculated next. At 
this point it is necessary to decide what part of the Z° —> ijj + £ + £~ phase space is to be 
identified as photon fragmentation, and what part is interpreted as lepton fragmentation. 

In Eq. (|9]), the dependence on the factorization scale \i cancels between and dT/dE^. 

Thus by changing fi some of the lepton fragmentation contribution can be moved into the 
photon fragmentation term and vice versa. There is therefore no clear crossover from lepton 
fragmentation to photon fragmentation, which makes it necessary to make an arbitrary 
choice on the appropriate phase space cutoff. In this paper a cutoff on the invariant mass 
of the £ — ip system, where £ is the fragmenting lepton, is introduced. The contribution to 
the differential decay rate from negative lepton fragmentation is considered to come from 
the region of phase space where s < fi 2 . Here s is the invariant mass of the £~ — ip system. 
Similarly the contribution from positive lepton fragmentation is considered to come from the 
region s' < fi 2 where s' is the invariant mass of the £ + — ip system. The photon fragmentation 
contribution is interpreted as coming from the remaining region. In the calculation of the 
photon subprocess the invariant-mass cutoffs translate into a limit on the phase space of the 
photon energy distribution. 
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The decay rate for Z° — > 7 + 



is 



T(Z° -> 7 + £ + r ) = ^ / [d*] [*+] [dp-] (2vr) V(Z — k — p + — p_) i |A| 2 , (18) 

where fe, p+, p_ and Z are the 4-momenta of the photon, £ + , and Z°. The amplitude 
can be calculated from the diagrams in figure 4. Averaging over initial spins and summing 
over final spins reduces the square of the amplitude to 



-Y\A\ 2 



9je 2 
6 cos 2 9„, 



,2 /( .(y-2) 2 + ™ 2 



(19) 



where y = 2k ■ Z/M 2 Zl and w = 2(p + - p_) • Z/Mf , C y = -1 + 4 sin 2 tt and C A = 1, ^ 
is the weak coupling constant, and 9 W is the weak mixing angle. Note that aside from the 
normalization this is the same as the integrand of Eq. (ffl) in the limit —>■ 0. Simplifying 
the phase space integral, Eq. Jl8| ) reduces to 



r(z°^ 7 + £ + r) 



M 9 . 



32(2vr) 3 hv*/M% 



dy 



(y-2^/M 2 z ) 



(20) 



where the limits on the invariant masses s and s' translate into the limits on w. Integrating 
over w gives the photon energy distribution in the decay Z° — > 7 + : 



(z°-. 7 (£ 7 ) + ^r,/i 2 ) 
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2a T(Z° 

7T 



M 9 



(y-l) 2 + l , ( v-i*IM% 
y H M 2 /Mi , 
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where y = 2E 1 /Mz- It is possible to simplify the expression for the energy distribution by 
taking the limit fi <C Mz- In this limit Eq. fl2~l"D simplifies to 



(Z°^7(^)+^-,/x 2 ) 



2« nz° -> £+r 



7T 



Ms 



(v 



l) 2 + 1 , yM\ 



log ; 



[22\ 



The price that is paid for this simplification is that smooth threshold behavior at y — 2/z 2 /Mf 
is lost, and the differential decay rate becomes negative at sufficiently small values of y. Since 
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the fragmentation approximation breaks down in the threshold region anyway nothing is lost 
by making this simplification. 

Lepton Fragmentation 

The calculation of the fragmentation function Dg^{z, /i 2 ) for a lepton to split into 
a ip parallels that of the photon fragmentation function. At leading order in a the process 
is symbolically represented by the diagram in figure 5a. The fragmentation probability is 
obtained by dividing o"i, the cross section for the production of ip + t~ , by 00, the cross 
section for lepton production shown in the diagram of figure 5b, in the limit 3> 
where fragmentation dominates. 

The general form of the cross section <To for lepton production is 

00 = -^j— J [dq][dp out ] (2ir) 4 5 4 {p in -q- p out ) ^ |A | 2 , (23) 

where q, Pi n , and p out are the 4-momenta of the l~ , the incoming particles, and all other 
outgoing particles. Just as in the photon fragmentation calculation, [dq] and [<ip «t] are the 
Lorentz invariant phase space for the lepton and the remaining outgoing particles, and Flux 
represents the incoming particle flux (which will cancel with the same quantity in the cross 
section o"i). The square of the amplitude A calculated from the Feynman diagram in figure 
4b for i~ production, averaged over initial spins and summed over final spins, is 

^|A | 2 = tr[^rf] (24) 

where the Dirac matrix T is the matrix element for the production of a real lepton of 
momentum q for which the explicit form is not needed. The lepton mass mi has been 
neglected since its 4-momentum q is taken to be large compared to mi. 

The general form of the cross section o\ for the production of a lepton that subse- 
quently fragments into a ip is 

— ^ J[dP][dpJ\[d Pout ] (2ir)U\p m -P-p^-p out ) EIAI 2 , (25) 

where P is the ip 4-momentum, and everything else is as described after Eq. (|23|). The next 
step is to write the phase space in an iterated form, by introducing integrals over q = P + p_ 

9 



the virtual lepton momentum, and s = q 2 the invariant mass of the if) — £ system. Then 
the phase space expression becomes 

[dP] [dp_] (27T) 4 5 4 ( Pm -P-p_- Pout ) 

= J^J \M (2vr) 4 5 4 fe„ - q - Pout ) J [dP}[dp_] (2n) 4 S\q - P - p.) . (26) 

The contribution that corresponds to the fragmentation of the lepton in the diagram of figure 
5a comes from the region of phase space in which the if) — 1~ system has large momentum q 
compared to the if) mass and small invariant mass s = q 2 of order M 2 . In a frame in which 
the virtual lepton has a 4-momentum q = (q , 0, 0, q^), the longitudinal momentum fraction 
of the if) relative to the ip — l~ system is z = (Pq + -P3) / {qo + qs) and its transverse momentum 
is P± = (Pi, P 2 ). Expressing the phase space in terms of these variables and integrating over 
the 4-momentum p_ and over P±, the 2-body phase space reduces to Jl 



2 - 



/] rl ( Ap \ 

[dP] [dp.] (2vr) V(g -P-p_) = -^dz6[s-^y (27) 

The lepton mass has been set to zero. An upper limit on the integral over the invariant mass 
s is introduced by requiring s < /i 2 , as discussed earlier. 

The calculation of the lepton framentation function is simplest to do in the axial 
gauge, because only the diagram of figure 5a, where the lepton is produced and splits into a 
collinear lepton and if), needs to be considered. Other diagrams where both the lepton and 
if) are produced separately from the vertex Y are suppressed in this gauge. If the calculation 
were done in some other gauge these diagrams would need to be considered, but the resulting 
expression could be manipulated into the form below using Ward identities. The 4- vector iV 
associated with axial gauge is chosen to be N = (1,0,0, —1). The amplitude A\ calculated 
in this gauge can be reduced to 

A, = e 2 ^ ep(Py - (u(p-) la ^T) [g^ ^ J . (28) 

where e is the if> polarization vector, and T is the Dirac matrix element for the production 
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of a virtual lepton with an invariant mass s on the order of the if) mass. The explicit form 
of T is not needed. 

The if) 4-momentum can be written as P M = zq^ + P± + (-Po<?3 — Q0P3) I {% + 93)-^- in 
the fragmentation region P± = (0, P±, 0), and Po<?3 — Q0P3 ar e of order while the virtual 
lepton momentum q is large compared to so P ~ zg. Note that in the fragmentation 
region both s = q 2 and P ■ q are of order M 2 . Using these approximations and keeping only 
the leading order terms in simplifies the square of the amplitude to 

ElA| 2 = 32 7 rV4( (z - 1 J + 1 l - ^) try TV]. (29) 

Integrating over s up to the scale fi 2 the lepton fragmentation probability is obtained by 
dividing the cross section o\ by the cross section do. The differences between <j\ and <To are 
due to the fact that q 2 ~ M 2 in a\, while q 2 = in a . These differences are on the order of 
M^/E^ and in the fragmentation limit where 3> M$ they can be neglected. The result 
is 

* = f is f iz - % 2o>4 f(»- 1 ) > + 1 I - (30) 
a JO Jo z v \ z s s 2 J 

from this it is possible to extract the lepton fragmentation function 



log —J - z + 



#(z/i 2 - M 2 ) (31) 



Note that the lepton fragmentation function is zero for values of z at which the production 
of a if) is kinematically forbidden. Taking the limit /i 3> simplifies Eq. ([Jl]) to 



D e ^(z, /1 2 ) « 2a 2 4 



l) 2 + 1 . z/i 



2 



(32) 



Just as in the calculation of the photon subprocess there is a price to be paid for this 
simplification. Eq. (^) does not have the correct threshold behavior at z = M 2 //i 2 , and 
it becomes negative for sufficiently small z. The fragmentation function Eq. (|32|) is shown 
at the scales /i = 3M^ and \x = 6M^ in figure 6. Note that there is a dramatic dependence 
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on the arbitrary factorization scale //. At the scale /i = 6M^, the fragmentation function is 
much larger and peaks at a lower z value than at the scale fi = 3M^. 



Comparison with full calculation 



The fragmentation contribution to the differential decay rate for Z' 



i> + £ + £~ is 



given by inserting P 7 _>^ from Eq. (|17D , dT/dE^ from Eq. (|22"D , and from Eq. (|32|) 



Note that the /z-dependence cancels exactly. It is now possible to verify that this agrees with 
the full calculation in the fragmentation limit Eq. (|j). 

Figure 7 compares the energy distribution of the full calculation Eq. (^) and the 
fragmentation calculation Eq. ([33]). It is clear from the graph that the fragmentation 
approximation breaks down for sufficiently small E^. For E^ = 3M^ the difference between 
the two curves is less than 1%, while for E^ = 2M^ the difference is 5%. In practice there 
will often be a minimum energy below which detectors do not register particles. If this 
minimum energy is large enough then it is clear that the fragmentation approximation will 
give a result very close to the full calculation. 

Figure 8 shows the energy distribution in the fragmentation limit separated into the 
lepton fragmentation contribution, the first term on the right hand side of Eq. (|^), and 
the photon fragmentation contribution, the last term on the right hand side of Eq. ([|). 
The contributions are shown at fi = 3M^, and fi = 6M^. The relative contribution of the 
two production mechanisms depends dramatically on the factorization scale n, though the 
total, photon fragmentation plus lepton fragmentation is independent of /i. At the scale 
fi = 3Mip the contribution of the lepton fragmentation mechanism is negligable compared 
to the contribution from the photon fragmentation mechanism, while at the scale [i = 6M<^ 
both contributions are of the same order. 

The analysis carried out in this paper applies equally to other heavy quark anti- 
quark states such as the ip' and the T. Unfortunately the framentation contribution to 



into the factorization formula Eq. (|D|): 



(Z° ^ij(E^)+ £+£-) = 8a 2 gl 



f(Z°^£+e-) [(y- 1)2 + 1 v *Mi 



2y . (33) 



dEf 
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T production is an order of magnitude smaller than the fragmentation contribution to ip 
production. This is because the T-photon coupling g\ = 6 x 1CT 4 is so much smaller than 
the ?/>-photon coupling — 8 x 10~ 3 , and the T mass is larger than the ip mass. 

Conclusion 

The process Z° — > ip + £ + £~ has been studied in the fragmentation limit Mz — > oo 
keeping E^/M z fixed. In this limit the decay rate factors into the subprocess rates T(Z° — > 
and f(Z° — > 7 + £ + £~), convoluted with the electromagnetic fragmentation functions 
and D 7 _»^. The fragmentation function D^^(z,/i) for a lepton to split into ip, the 
fragmentation probability P 7 _>^, for a photon to split into if), and the subprocess T(Z° — > 
7+£ + £^) were calculated at lowest order in a. The lepton fragmentation function was defined 
by imposing a cutoff on the invariant mass s of the lepton and tp in the final state. This cutoff 
translated into limits on the phase space of the subprocess T(Z° — * 7 + £ + £~). It was then 
explicitly shown that the /x-dependence of the lepton fragmentation function canceled the 
/^-dependence of the subprocess T(Z° — > 7 + £ + £~). Comparison between the fragmentation 
calculation and the full calculation shows that the fragmentation approximation is accurate 
to within 5% at = 2M^, and becomes more accurate for ip energies greater than this. 
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Energy Physics, under Grant DE-FG02-91-ER40684. I wish to thank E. Braaten for many 
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Figure Captions 

1. The Feynman diagrams for Z° — > vp + £ + £~ at leading order in a. a) The two diagrams 
from which the fragmentation contribution can be isolated, and b) the two diagrams 
that may be neglected. 

2. The two Feynman diagrams for Z° — > ip^y at leading order in a. 
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3. Feynaman diagrams for a) ip production by photon fragmentation, b) photon produc- 
tion. The shaded circle represents some general vertex that radiates the photon. 

4. The two Feynman diagrams for Z° — > 7 + £ + £~ at leading order in a. 

5. Feynman diagrams for a) ip production by lepton fragmentation, b) lepton production. 
The shaded circle represents some general vertex that radiates the lepton. 

6. Lepton fragmentation function for /x = 3M^ (solid) and /1 = 6M^, (dashes). 

7. Energy distribution: the full calculation (solid), the fragmentation calculation (dashes). 

8. Photon fragmentation (PF) and lepton fragmentation (LF) contributions to the energy 
distribution for /j, = 3M^ (solid) and (j, = QM^ (dashes). 
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